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Abstract 

Paramagnetic effects are shown to result in the appearance of a triplet component of order 
parameter in a vortex phase of a d-wave superconductor in the absence of impurities. This 
component, which breaks both parity and time-reversal symmetries of Cooper pairs, is expected 
to be of the order of unity in a number of modern superconductors such as organic, high-T c , 
and some others. A generic phase diagram of such type-IV superconductors, which are singlet 
ones at H = and characterized by singlet-triplet mixed Copper pairs with broken time-reversal 
symmetry in a vortex phase, is discussed. 
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It is well known [1,2] that type-II superconductors, where superconductivity survives at 
high magnetic fields, H cl (T) < H < H c2 (T), as Abrikosov vortex phase [3,2,1], are subdi- 
vided into two main classes. They are superconducting alloys (or dirty superconductors) 
[1,2] and relatively clean materials, where type-II superconductivity is due to anisotropy of 
their quasi-particles spectra and relatively heavy masses of quasi-particles [4]. The latter 
compounds are currently the most interesting and perspective superconducting materials, 
including organic [5], heavy fermions [6], high-T c [7], MgB 2 [8], and some other supercon- 
ductors. 

Superconducting orbital order parameter, A(r 1; r 2 ), corresponding to pairing of two elec- 
trons in Cooper pair, can usually be expressed as A(ri,r 2 ) = A(R)A(r) [9,10], where 
external order parameter, A(R), is related to motion of a center of mass of Cooper pair, 
R = (ri + r 2 )/2, whereas internal order parameter, A(r), describes relative motion of 
electrons in Cooper pair, r = r x — r 2 . In this context, type-II superconductors in their vor- 
tex phases are characterized by broken symmetries of A(R), corresponding to vortices and 
Meissner currents. Other important issues are symmetries of internal orbital order param- 
eter, A(r), and related spin part of order parameter, A(a 1 ,a 2 ). To satisfy Fermi statistics, 
in case of singlet superconductivity (where the total spin of Cooper pair |S| = 0), internal 
order parameter, A(r), has to be an even function of coordinate r, whereas, in case of triplet 
superconductivity (where |S| = 1), A(r) has to be an odd function of r. In accordance 
with symmetry properties of A(r) (or its Fourier component, A(k)), superconductors are 
subdivided into conventional ones [1,2] (where superconductivity is described by BCS s-wave 
singlet pairing) and unconventional ones [9,10] (where symmetry of A(k) is lower than un- 
derlying symmetry of crystalline lattice). At present time, unconventional d-wave singlet 
superconductivity is firmly established in high-T c [11] and some organic materials. There 
exist also several strong candidates for unconventional triplet superconducting pairing such 
as Sr 2 Ru0 4 [12], (TMTSF) 2 X [13], ferromagnetic [14], and heavy fermion [9,10,15] supercon- 
ductors. It is a common belief that magnetic field does not change internal order parameters, 
A(k) and A(<7i, <r 2 ), in conventional [1,2] and unconventional [9,10] type-II superconductors. 
Moreover, although Meissner currents break time- reversal symmetry of A(i2), internal order 
parameters, A(k) and A(<ti,<t 2 ), are believed to preserve t — > —t symmetry. 

The goal of our Letter is to demonstrate that there have to exist type-IV superconductors 
[16], which exhibit qualitatively different magnetic properties. More precisely, we suggest 
and prove the following theorem: each singlet type-II superconductor in the absence of 
impurities is actually type-IV superconductor with broken k — > — k and t — > —t symmetries 
of internal Cooper pairs wave functions in vortex phase, provided that effective constant of 
triplet pairing is not exactly zero, g t ^ 0. We show that the above mentioned theorem is an 
inherent property of singlet superconductivity and is due to careful account for paramagnetic 
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spin-splitting effects in vortex phase, which have been treated so far only for g t — [17,1]. 

We define type-IV superconductivity as singlet superconductivity at H = which exhibits 
broken symmetries of internal Cooper pairs wave functions, A(k) and A(<7i,<72), in vortex 
phase. In our particular case, internal order parameter in vortex phase is shown to be a 
mixture of a singlet d-wave component, A s (k), with a triplet component, i A t (k). Note that 
this order parameter breaks not only parity symmetry, k — > — k, but also a time-reversal 
symmetry, t — > —t, due to an imaginary coefficient i. Below, we demonstare that the effects 
of singlet-triplet coexistence are of the order of unity in a number of modern clean type-II 
superconductors, where the orbital upper critical fields are of the order of paramagnetic 
limiting fields, fisH^O) ~ T c [18] (see Table 1). It is important that the suggested theorem 
is very general: it is valid even for simplest spin independent electron-electron interactions 
for both attractive and repulsive interactions in triplet channel. As discussed below, this 
theorem is based on symmetry arguments and is a consequence of broken spin symmetry (due 
to paramagnetic effects) and broken translational invariancy of A(R) (due to the existence 
of vortices). 

Here, we discuss how paramagnetic effects result in the appearance of a triplet component 
in vortex phase using qualitative arguments. We recall that, in vortex phase, external order 
parameter, A(R), is a function of R on scale of £, where £ is a coherence length [1-4]. There- 
fore, A(R) corresponds to superconducting pairing of electrons with total non-zero momenta 
of Cooper pairs of the order of |q| ~ %/£,. As seen from Fig.l, a probability of pairing for 
electrons with spin up and spin down, |A(+, — )| 2 , is different from that for electrons with 
spin down and spin up, |A(— ,+)| 2 , if q 7^ 0. Therefore, singlet superconductivity, which is 
characterized by spin order parameter A(+, — ) = — A(— , +) [9,10], has to be mixed with a 
triplet component, characterized by spin order parameter A(+, — ) = A(— , +) (i.e., |S| = 1 
and S y = 0) (see Fig.l). 

Below, we quantitatively describe superconducting pairing with internal order parameter, 
exhibiting broken inversion and time-reversal symmetries, in d-wave singlet superconductor 
with layered electron spectrum, 

e (k) = (kl + k 2 y )/2m + 2t± cos(M) , e F = mv 2 F /2 , (1) 

in a magnetic field: 

H = (0,#,0), A = (0, 0, —Hx) . (2) 

In case, where electron-electron interactions do not depend on electron spins, the total 
Hamiltonian of electron system can be written in the form: 

H = H + H int , Ho = E e -( k ) a * ( k ) a -( k ) > 

k,<7 



3 



Hint = ^EE VQl, kO a+(k + S) a + CT (_k + 5) a _ ff (_ kl + 2) a CT ( kl + 5) , (3) 

q,<r k,ki 

where e CT (k) = eo(k) — (jubH (a = ±1), a+(k) and a ff (k) are electron creation and annihi- 
lation operators. As usually [9,10], electron-electron interactions are subdivided into singlet 
and triplet ones: 

V(k, k x ) = V s (k, kx) + H(k, ki) , V s (k, kx) = K(-k, kx) = V s (k, -ki) , 

\4(k, ki) = -V t (-k, ki) = -V t (k, -ki) . (4) 

We define normal and Gorkov (anomalous) finite temperature Green functions, 

G a , a (k, ki; r) = - < T r a CT (k, r)a+(ki, >, F a _ CT (k, k x ; r) =< T r a CT (k, r)a_ CT (-ki, 0) > , 
F+_ CT (k, ki; r) =< T r a+(-k, r)a^(ki, 0) > , (5) 

as well as singlet and triplet superconducting order parameters, 

A s (k, q) = ~ E K(k ki)T£[F + ,_(^ n ; kj + |, ki - S) _ F_, + (^ n ; k x + |, k x - |)] , 

A t (k, q) = ~ E ki)rE^,-K; k! + |, ki - |) + F_,+ (^ n ; k x + |, k x - |){Q) 

ki 

by standard ways [19,20,9,10]. 

The goal of our Letter is to consider a phase transition line between metallic and singlet- 
triplet mixed superconducting phases in Ginzburg-Landau (GL) region, (T c — T)/T c <C 1 
[1-4], where T c is a transition temperature between metallic state and d-wave singlet phase 
at H = 0. For this purpose, we linearize Gorkov equations [9,10,19] with respect to order 
parameters (6) and obtain the following system of linear equations [21]: 

A s (k, q) = ~ E ^( k > ki)T]T[A s (ki, q) S + A t (k 1; q) D] , 

kl <^>n 

A t (k, q) = ~ E ^( k ' ki)T E[A*(ki, q) S + A s (k x , q) £>] , 

ki w n 

5 = G° + (iu n , k x + |)G°_(-ia; n , -k x + |) + G°_(^ n , k x + (-ic; n , -k x + |) , 
L> = G° (iw n , k x + |)G°_ (-iw n , -kx + |) - G°_(iu; n , k, + ^)G° + (-icu n , -k, + |) , (7) 

where G®(iu n ,k) = l/(iu n — e CT (k)) is Green function of a free electron in the presence 
of paramagnetic effects and oo n is Matsubara frequency [20]. [Note that common Eqs.(7) 
directly demonstrate singlet-triplet coexistence effects in vortex phase, where q ^ (see 
Fig.l)]. 
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Below, we consider in detail an important example, coexistence of singlet d x 2_ y 2-wave [21] 
and triplet p^-wave order parameters, which corresponds to the following matrix elements 
of electron-electron interactions: 



K(k,ki) 
V t (k,kO 



47r / g s cos 2(f) cos 20! 



vf 



g s > 0, g s > g t , 



(8) 



g t cos(0-0i) 

where and (f>\ are polar angles corresponding to momenta k and ki, respectively. [Note 
that inequalities g s > and g s > gt guarantee that singlet d x 2_ y 2 -wave phase is a ground 
state at H — and T < T c ]. After substitution of Eqs.(8) in Eqs.(7), we represent order 
parameters as follows, A s (k, q) = v^2cos20 A s (q) and A t (k, q) = a/2cos0 A t (q), and 
rewrite Eqs.(7) in a matrix form: 



A ss (q) A st (q) \ / A s (q) \ = / A s (q)/g s 
A ts ( q ) A tt (q) ) { A t (q) J { A t (q)/gt 



(9) 



We calculate matrix A(q) at g y = in GL region [3,4,22,9] which corresponds to its 
expansion as power series in small parameters, v F q x jT c <C 1 and t±dq z /T c < 1. As a result, 
we obtain: 

/ 



A 



(2vrT) EL„>o 



-H B Hv F q x {>KT c ) J2Z>o 7? 



V 



-H B Hv F q x (TrT c )T,Z>o^ , 



(27TT) EC n >o 



i- ^(3^/2 + 4^) 

(10) 



with f2 being a cut-off energy. Magnetic field (2) is introduced in Eqs.(9),(10) by means of a 
standard quasi-classical approximation [23,22,4,3], q x — > —i(d/dx), q z /2 — > eA z /c = eHx/c, 
which leads to the following matrix GL equations extended to the case of triplet-singlet 
coexistence: 



I dx 2 



V 



* V2-K 



(2^x) 2 



2^2 



r + 



gt-g s _|_ 3^2 ( _d? 



\ 



J 




0, (11) 



where r = (T c - T)/T c < 1, = ^7C(S)v f /AV2ttT c and £ ± = ^/7C(3)t ± rf/2 v / 27rT c are 
GL coherence lengths [4,9], £(3) ~ 1.2 is zeta Riemann function, 0o is a flux quantum, x is 
coordinate of a center of mass of Cooper pair. In typical case, where g s — g t ~ <? s , Eqs. (11) 
have the following solutions: 



A s (x) 

Mx) 



exp - 



2 « 



V7C(3) 



exp 



TX 



(12) 



Eqs.(ll),(12) are the main results of our Letter. They extend GL differential equation 
[1-4,22,9,10] and its Abrikosov solution for superconducting nucleus, exp(— rx 2 /2£j) [3,1,2], 
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to the case g t ^ 0. Eqs.(ll),(12) directly demonstrate that, in vortex phase, singlet order 
parameter always coexists with triplet one, characterized by |S| = 1 and S y — (H | 
y), for arbitrary sign of effective triplet coupling constant g t . It is important that triplet 
component (12) breaks not only parity but also time-reversal symmetry due to the existence 
of non-diagonal matrix elements, proportional to iH, in Eqs.(ll). Indeed, A£(x) ^ A t (x) in 
Eq.(12), which indicates that t — > —t symmetry is broken and, thus, Cooper pairs possess 
some internal magnetic moments [9,10]. 

To summarize, the main message of the Letter is that Cooper pairs cannot be considered 
as unchanged elementary particles in vortex phases of modern strongly correlated type-II 
superconductors, where fiBH c2 (0) ~ T c and \g s \ ~ \g t \. Indeed triplet-singlet components 
ratio in Eqs. (12) at x — y/2£\\/y/r and low temperatures, r ~ 1, can be estimated as 
R = A t /A s ~ % [fi B H C 2(0)/T c ] (see Table 1). The appearance of a triplet component, 
breaking time-reversal symmetry, has to change all qualitative features of vortex phases in d- 
wave superconductors. These include the appearance of non-zero internal magnetic moments 
of Cooper pairs, possible unusual topology of superconducting vortices, the appearance of 
spin-wave excitations, the disappearance of quasi-particles near zeros of d x 2_ y 2-wave gap, 
cos 20 = 0, unusual spin susceptibility, and other non-trivial phenomena to be studied in 
the future. We suggest that, in clean type-II superconductors, there exist the forth critical 
fields, H ci (T), corresponding to phase transitions (or crossovers) between Abrikosov vortex 
phases and some exotic vortex phases and call such materials type-IV superconductors. In 
conclusion, we point out that singlet-triplet mixing effects were earlier studied in He 3 [24], 
Larkin-Ovchinnikov-Fulde-Ferrell phase [25,26], for surface superconductivity [27], and in 
superconductors without inversion symmetry [9,28,29]. 

The author is thankful to N.N. Lebed (Bagmet), P.M. Chaikin, and V.M. Pudalov for 
numerous and useful discussions. 
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TABLE I: Upper critical fields, H c2 (0) [30-33], transitions temperatures, T c , and triplet-singlet 
ratio, R, are listed for some modern layered d-wave and s-wave superconductors. 



- {ET) 2 AuI 2 P - {ET) 2 IBr 2 YBa 2 Cu 3 7 MgB 2 

H c2 (0) [T] ^OD 2^(jjj H0(I) 18(H) 

VbH c2 (0) [K] 3.7 1.6 74 12 

T c [K] 4.3 2.3 85 35 

R 0.85 0.7 0.85 0.4 
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FIG. 1: Paramagnetic effects split electron spectra with spin up and spin down: e + (k) = eo(k) — 
HbH and e~(k) = eo(k) + fJ-sH, correspondingly. Two Cooper pairs with spin parts of wave 
functions, A(+, — ) and A(— ,+), and equal total momenta, q 7^ 0, are characterized by different 
probabilities to exist since energy difference — e±\ = qvF + ^I^bH is not equal to energy 
difference — = —qvF + ^HbH . [For simplicity, linearized one-dimensional electron spectrum, 
e(k) = VF\k\, is shown]. 
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